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The quantum processes involved in the interaction of matter with, separately, an electron vortex 
(EV) and an optical vortex (OV) are described, with matter modelled in terms of a neutral two- 
particle atomic system, allowing for both the internal (electronic-type) motion and the gross (center 
of mass- type) motion of matter to be taken into account. The coupling of the atomic system to 
the EV is dominated by Coulomb forces, while that of the OV is taken in the p • A canonical form 
which couples A, the transverse vector potential of the optical vortex, to the linear momenta of 
the two-particle system. An analysis of the dipole active transition matrix element is carried out in 
each case. The electron vortex is found to be capable of exchanging its orbital angular momentum 
(OAM) with both the electronic and the center of mass motions of the atomic system in an electric 
dipole transition. In contrast, for electric dipole transitions the optical vortex is found to be capable 
of exchanging OAM only with the center of mass. The predictions are discussed with reference to 
recent experiments on electron energy loss spectroscopy (EELS) using EVs traversing magnetised 
iron thin film samples and those involving OVs interacting with chiral molecules. 
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I. INTRODUCTION 



Vortices are now well established phenomena in a number of physical contexts, most notably in condensed matter 
physics, where they are known to play an important role in the dynamics of superfluid flow [l| and, more recently, 
in dilute gas Bose-Einstein condensates 0; in optics as optical vortices in the form of Laguerre-Gaussian and Bessel 
light beams 041] and, very recently, in electron microscopy as electron vortices Broadly, a vortex field has 

a propagating wavefront endowed with a screw dislocation. The distinctive feature throughout is the presence of a 
phase factor exp (il(f>) where <f> is the azimuthal angle about the beam axis, and I, an integer taking positive and 
negative values, is the winding number, such that the field posesses orbital angular momentum (OAM) of ITi. Over 
the last two decades or so, much work has been carried out on optical vortex (OV) beams, so much so that optical 
vortex physics has now been established as a new branch of modern optics. The area began with the seminal work 
by Nye and Berry [lj| who put forward the suggestion of beams endowed with the vortex property, and the work by 
Allen et al. [ll[, first reported in 1992, which set the scene for the subsequent research into optical vortices and their 
interaction with matter. Optical vortex physics is rich in both fundamental considerations and applications, and has 
led to key applications in the manipulation of matter, both in the bulk and at the level of basic constituents. This 
impact is envisaged to continue in inter-disciplinary areas, including quantum information processing. 

The phenomenon of electron vortices (EVs) is a very recent addition to the growing catalogue of vortex physics, 
essentially arising from the concept of OVs. The suggestion for their existence was first made by Bliokh et al. [l|. 
Like an OV, an EV also carries the key property of quantized OAM of Yh per electron. Following Bliokh et al.'s 
suggestion, experimental work led to the creation of EVs, first by Uchida and Tonomora 0] using a stepped spiral 
phase plate, followed by Verbeeck et al. Q who used a binary holographic grating with a Y-like point defect, and 
thirdly by McMorran et al. [t| who also used the holographic plate technique but managed to generate EV beams 
with winding numbers as high as I = 100. It is now clear that EVs can be readily generated inside an electron 
microscope. More recent advances include the creation of EV beams of cross sections in the atomic scale [H[ and 
the realisation that they can also be generated using spiral phase plates [l4[. EV beams arc predicted to lead to 
important new advances in the physics and potential inter-disciplinary applications of electron beams, revolutionizing 
electron microscopy and spectroscopy with atomic scale resolution in the imaging of materials, including those with 
low absorption contrast such as biological specimens (l5j . 

Since OAM is a well defined quantized property both of OVs and EVs, it is natural to expect it be exchanged when 
cither type of vortex beam interacts with matter in the form of atoms, molecules or solids. Both optical spectroscopy 
using OV beams and electron energy loss spectroscopy (EELS) using EV beams are expected to involve an exchange 
of quantized OAM with matter, especially in the normally dominant electric dipole transitions. The exchange must 
involve both the 'electronic-type' motion and the 'center of mass-type' motion of matter. The purpose of this paper 
is to systematize the theory appropriate for the processes of OAM exchange in the interactions of OVs and EVs, 
separately, with a model system of matter in the form of a neutral two-particle system. The main aim is to find out 
whether, and if so, in what manner, do EVs differ from OVs in the processes of exchange of OAM. This is done by 
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extracting information residing in the respective transition matrix elements leading to the selection rules associated 
with OAM transfer. 

This paper is organised as follows. In section |H] we outline the basic formalisms for the OVs and EVs in the forms 
of free Bessel modes, define the model atomic system and emphasise the necessity of separating the matter motion 
into internal (electronic-type) and gross (center of mass-type) motions. We write down the total Lagrangian of the 
system relative to the laboratory frame and follow canonical steps to determine the conjugate momenta, leading 
to the total Hamiltonian including the respective interactions, which we then restrict to the dipole approximation. 
Finally we determine the unperturbed quantum states participating in a typical transition. In section Hill we consider 
OVs in interaction with the model atomic system and proceed to discuss the selection rules governing the dipole 
transition matrix element. We also discuss the experimental results emerging from work by Araoka et al. [l6[ which 
conforms with the theoretical prediction for this case. In section [TVl we consider the corresponding theory for EVs 
interacting with the model atomic system and deduce the selection rules governing the dipole transitions in this case 
as well. In section [V] we discuss the predictions of our theory in relation to the results of the recent experiment 
on EELS by Verbeeck et al. [|[ who reported the observation of dichroism in electron energy loss spectroscopy of 
thin film magnetised iron samples using EVs. Section I VII contains a summary of the main conclusions regarding 
the comparison between the two vortex interactions with matter and provides further comments. Some details are 
consigned to Appendices. 



II. OPTICAL AND ELECTRON VORTEX MODE FUNCTIONS 



A. Bessel modes 



For both types of vortices we shall concentrate on the simplest type of vortex function, namely the Bessel mode. As 
a solution to both the Hclmholtz and Schrodinger equations the Bessel beam is a suitable carrier for orbital angular 
momentum for both optical and electronic systems; choosing solutions involving only Bessel functions of the first kind 
gives a mode with zero intensity along the optical axis, as required to support the phase singularity arising from the 
vortex phase factor e tl ^. 

In the first case of the optical vortex the mode function is characterised by a transverse electric field which is a 
solution of the electromagnetic vector Hclmholtz equation. In cylindrical polar coordinates r = (p, </>, z) the optical 
vortex mode function is 

E(r,f) - E J l (k ±P )e tk ' z e u ' p e- luJt i, (1) 

where Eq is the mode amplitude and the unit vector e denotes the wave polarisation vector. Since we are mainly 
concerned with orbital angular momentum, we shall assume that e stands for linear polarisation, unless stated 
otherwise. The radial function Ji(k±p) is the Bessel function of the first kind of order I, where I is the winding number. 
The wavevectors fcj_ and k z stand for in-plane and axial wavevector variables respectively, such that k 2 = k\ + k 2 , 
and ui is the frequency of the light. For the optical case, the transverse vector potential is related to the vortex electric 
field of Eq. Q by 

A(r ) t) = --E(r,t) (2) 

The electron vortex mode is characterised by a wavefunction ip(r,t), which is a solution of the scalar Hclmholtz 
equation emerging from the Schrodinger equation in cylindrical polar coordinates. For convenience and for ease of 
comparison we retain the same symbols, including wavevector and frequency variables. We have 

VM) = JVJ|(fcip)e lt2 e ,! *e- rat , (3) 

where N is a suitable normalisation constant and cj = £/h with £ the mode energy. For both types of vortex 
mode the vorticity resides in the phase factor e l1 ^. Note that the two vortex functions described above are similar 
in appearance and formally have the same spatial and temporal distributions; however their physical characteristics 
differ markedly, firstly in the scales of variation and, secondly, since they describe very different phenomena, their 
respective interactions with matter differ significantly. For the optical vortex case, the coupling is via the minimal 
coupling prescription, leading to interaction terms of the form p • A(r, t) representing the interaction of the transverse 
vortex vector potential with the momentum, p, of each of the atomic particles. In contrast, the leading interaction of 
the charged electron vortex is a coupling to the two-particle atomic system via the Coulomb interaction, with scalar 
potential $(r') coupling to each of the atomic constituents. We have 

$(r') = (4) 

* [r> 47re K -r'| { ) 
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FIG. 1. The relevant coordinate frames in the description of the interaction between a two-particle neutral system and Bessel- 
type optical or electron vortex beam (schematic). The vortex position variable, r„ relative to the laboratory frame is given 
in cylindrical coordinates ; R is the position variable of the atomic center of mass, and q stands for the position variable of 
the internal (electron-type) motion. The projections of the three position vector variables on the xy plane are seen to have 
azimuthal anglular positions 4> v , cj>R, and cj> q respectively. 



as the Coulomb potential at r' due to the electron vortex with position variable r t , . 

The main task is to explore, by direct analysis, in each case how the transfer of OAM can occur in transitions between 
states of the vortex beam and the two-particle system, and deduce the selection rules governing those processes. To 
this end we use a rigorous canonical theory based on the Lagrangian which leads us to the total Hamiltonian of the 
unperturbed states and their interactions. 



B. Total Lagrangian and Hamiltonian 



With the laboratory coordinate system taken to be centered at the vortex origin, as shown in Fig. [T] we shall 
assume that the two-particle system is hydrogenic, consisting of an electron of charge — e, mass m e , position vector 
r e and momentum p e ; and a nucleus of mass m p , charge e, position vector r p and momentum p p - all relative to the 
(laboratory) coordinate system in which the vortices are defined with position variable r„ , and such that all position 
vectors are in cylindrical coordinates r a = (p a , (j) a ,z a ). In this frame of reference the two-particle system possesses a 
total charge density given by 



p(r) = eS(r - r p ) - eS(r - r e ) 



(5) 
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The Lagrangian of the two-particle system in the presence of the interaction with both vortex beams can be written 
as 

L = -^atomic + ^vortex + ^ (6) 

Explicitly we have (following standard Lagrangian techniques as in, for example (l7j ) 

1 1 e 2 I 
^atomic = oW e re + + 7ZTTZ Z~i ' ( 7 ) 



-.tame - 2 »^e T- 2 »^p-r ^ 



Cortex = f / (A 2 (r,t)-c 2 (VxA(r,i)) 2 )A + ^,; (8) 

L int = er p ■ A(r„) - er e ■ A(r e ) - J p(r')$(r')d 3 r' . (9) 

Thus ^atomic i s simply the sum of the kinetic energies minus the inter-particle Coulomb potential energy of the atomic 
particles, and L° ortox is the zero order energies of the optical and electron vortices. L lnt contains the interaction of 
the optical vortex vector potential with each of the atomic particles, as well as the interaction of the electron vortex 
Coulomb potential with the charge density of the atomic particles. 

The formalism so far is such that the coordinates of the two-particlcs forming the atom are expressed entirely 
relative to the laboratory frame. However, to be able to discuss transitions involving the internal states of the atom, 
we proceed by expressing the atomic motion in terms of the gross (center of mass-type) motion and the internal 
(electronic-type) motion. This is achieved by the following transformation 

q = r e -r p ; R= M P P , (10) 

where M = m e + m p is the total mass. We now have q as the internal coordinate representing electronic-type motion 
about the nuclear position r p , while R is the coordinate of the center of mass in the laboratory frame. In terms of q 
and R the particle position variables are given by 

r e = R+-^q; r p = R - ^q (11) 

Substituting for r e and r p using Eq. (fTTj) . the Lagrangian becomes 

1 / • m. D . \ 2 1 / • m e . \ 2 / • m v . \ . , . / • m 



i= 2 TOe (^ + ^) + 2 W ^-^) " e (*U ^4) -AW + e (ft -^4). A(r,) 



,2 i „2 



1 e z 1 



4ire Q q| 47re 



-^vortex (12) 



where, for convenience at this stage, we have left unchanged the implicit dependence of the vector potential A(r) on 
the original particle coordinates, and have evaluated the integral involving $(r'). We note again that r v denotes the 
vortex position variable. We now determine the momenta pr and p q canonically conjugate to R and q respectively. 
We have after rearrangements 

dL 

PR = ^ = MA-e(SA) (13) 



dL 

dq 

where we have introduced the short hand notation 



P< Z = ^T = Mq-e(SA), (14) 



(SA) = A(r e ) - A(r p ) (15) 
(SA)^A(r e ) + ^A(r p ). (16) 

The momentum canonically conjugate to is p„ = m e r v and the momentum density canincally conjugate to A is 
II = eoA. The Hamiltonian representing the two-particle atom in the presence of the interaction with the two vortices 
is given by the standard expression 

H = p R R + • q + p v -r v + / f n ■ A) d 3 r - L. (17) 
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Substituting for all the canonical momenta and momentm density we find that a number of terms cancel to yield the 
result 

H = -MIL 2 + l M q 2 ( - -) + iJ° (18) 

2 2^* 47re |q| 47re V K - r p | [r^ — r e | ^ v ' V ' 

where fi = m e m p /M is the reduced mass and H® is given below in Eq. (f2"T)l . That the Hamiltonian in Eq. (fT5)) 
incorporates the interactions with the two vortices can be seen on eliminating the particle velocities in favour of the 
canonical momenta using. We have 

[p, + c(l]A)] 2 [p fl + e(£A)] 2 _ _ _e>_ ( 1 1_\ 

2 M + 2M 47reo |q| 4^0^, -r p | |r.„-r e |; + v " 1 9J 

Expanding the squares in the first two terms, we find 

H = ^ ~—-L+vk + £ p (SA) + — p R ■ (SA) -—( , r - , r^l + H° (20) 

2a* 47re |q| 2M /i Pg V ' M VR K ' ^e \\r v -t p \ |r„ -r e |y v V ' 

where we have retained only the terms linear in A(r), since for the optical vortex case our main concern is with 
transitions involving single optical vortex photons. Finally, we concentrate on the electric dipole approximation and 
expand the functions of (R — ^fq) and (R + -jf q) in Eq. (|20[) about the center of mass coordinate R, to find 

1 1 _ q ■ (r„ - R) , ^ 2 



|r,, -rJ |r,, -r R | |r„ - R| 3 



<W); (21) 



A(r e ) « A(R) + ^(q ■ V)A(R) + 0(q 2 ); (22) 

A(r p ) « A(R) - ^(q ■ V)A(R) + 0(q 2 ), (23) 

Retaining only those terms up to first order in q, as is consistent with the dipole approximation, we can then write 

H = H° q + H° cm + + + + (24) 

where the individual terms are expressed as 

o P 2 e 2 1 

H « = i~^H' (25) 

HL = (26) 
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K = jJ [A'(t) + c 2 (V x A(r))^ d A r + (27) 

^ov 9) = -P 9 -A(R); (28) 
A ^ 

tfov (i?) = ^P*-(q-V)A(R); (29) 



int _ e 2 q ■ (r„ - R) 
^ EV "4^ |r-R|3 (30) 

to dipole order. We can now write the overall Hamiltonian of the system as the sum 

H = H Q + H^ + H£l (31) 

where H° = H®+H® m +H® is the zero order Hamiltonian of the overall system , while Hq V and H^ v are the interaction 
Hamiltonians of the two-particle system with the optical vortex (including the Rontgen interaction given by Eq. (1291) , 
c.f. 18[) and the electron vortex, respectively. The total zero order Hamiltonian H° consists of four separate zero order 
Hamiltonians representing the distinct sub-systems, namely, H® representing the internal (electronic-type) motion, 
H® m representing the gross (center of mass-type) motion while, as mentioned above, H® = Hq V + H^, v represents the 
zero-order Hamiltonian of the free optical vortex and the electron vortex, respectively. We have 

H = H q + H cm + Fqv + h ev- (32) 
In this paper we will be dealing with only one vortex at a time interacting with the model atomic system, so that in 

'EV 



applications of Eq. (|32[) only three zero-order Hamiltonians are required. For the optical vortex case -H^y is excluded, 



while in the electron vortex case -Hqv * s excluded. 
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C. Unperturbed quantum states 

The selection rules can be deduced from a careful analysis of the transition matrix element evaluated between 
appropriate final and initial states [tyf } and \^> l ) of the vortex plus the two-particle system. These are eigenstates of 
the unperturbed Hamiltonian H°. For the optical case the initial and final states are products of the atomic electron 
and center of mass eigenstates, and the number state of the optical vortex field: 

|^ / > = |^ ,/ (q);^ / ( R );» <,/ ). ( 33 ) 

The number states ^'J) are eigenstates of the optical vortex Hamiltonian for which the vector potential operator is 
written in terms of annihilation and creation operators, ak ± ,k z an d fc as follows 

A(r,t) = A(r, t)a k± , kz + A*(r, t)al ±>kz . (34) 

For the electron vortex case, the initial and final states are written as product states of the electron vortex wavcfunction 
and the atomic internal electron and center of mass eigenstates, given as 

|^/) = |^/(q);^/(R) ; ^/(r,)). (35) 

The internal states of the atom \%/j q (q,9 q ,4>q)) = \n q ;£;m) are eigenstates of H® and can be formally identified as 
the well known hydrogenic states, here given in spherical coordinates (q,6 q ,(j) q ). Explicitly we have 

|^(q)) = \M9,0 q ,M) = N nAm Q n (q)PP(co S (e q ))e im ^ (36) 

where the integer I is the internal atomic orbital angular momentum (not to be confused with Z, the vortex OAM 
quantum number about the beam axis); m is the internal atomic magnetic quantum number (such that —£ < m < £), 
and n is the principal quantum number of the internal atomic motion. 

The eigenstates of the center of mass are taken to be product states of both its translational and rotational motion 

|Vfl(R)> = \MprAr,zr)) = H( PR )e' K ^e* K ^e lL ^. (37) 



where the subscript R indicates center of mass coordinates relative to the laboratory frame. Kr and K z are center 
of mass wavevectors for the in-plane translational motion and motion along the z-axis, such that the total linear 
momentum of the center of mass is given by K 2 = K 2 + K R . L is the orbital angular momentum quantum number 
of the center of mass about the beam axis. 

The vortex states of both the optical vortex and the electron vortex are described above by Eq. ([T]) and Eq. ((3]) 
respectively; vectors and their components relating to the vortex beams will henceforth be denoted with the subscript 
v. In Eq. (|33p and Eq. (|35[) , and also in what follows, initial and final values of functions or constants are respectively 
denoted by the superscripts i or /, while final values of quantum numbers are indicated by the presence of a dash. 

We seek to determine the orbital angular momentum selection rules in processes involving dipole active transitions 
due to the interactions between the vortices and the atom such that an exchange of orbital angular momentum occurs 
between the three subsystems. This will be done separately for the optical vortex and the electron vortex, in each case 
by analysing the complex matrix element (^|i?Q* E . )V |\E' 1 ) and its modulus square; the latter enters the well known 
formula for Fermi's golden rule leading to the evaluation of the transition rate. Note, however, that the orbital angular 
momentum is described by the azimuthal features of the optical vortex field and electron vortex wavefunction. For 
this reason only these azimuthal angular components in the matrix elements need be explicitly evaluated. In this way, 
the orbital angular momentum selection rules will be made apparent, and the question of transfer of orbital angular 
momentum between the vortex and the atomic internal dynamics will be answered. 



III. THE OPTICAL VORTEX CASE 



The two terms identified as dipole order interaction Hamiltonians, for the OV are given by Eq. ([25]) and Eq. (|29p . 
It can be shown that iJgy R ^ is in fact the same as the so-called Rontgen interaction Hamiltonian [l8| which couples 
the center of mass motion and the internal motion via the optical vortex magnetic field B. We shall not consider 
the effects of this Rontgen interaction Hamiltonian any further here, as it is typically much smaller than the direct 
coupling with the vector potential , given by Hq^^ (l9| . We now write Hq^^ in operator form: 

£$«>=-£p,.A(R). (38) 



7 



The interaction is proportional to the linear momentum operator p q of the internal atomic motion, and it should be 
noted that the vector potential operator is evaluated at the center of mass coordinate R, expressed in the laboratory 
frame of reference (see Fig. [I}. The transition matrix element is given by 



M 



fi 

ov 



^(q);^(R) 



H, 



int(g) 
OV 



^(q);^(R);; 



(39) 



That the interaction Hamiltonian contributes to the electric dipole transition by virtue of the presence of internal 
momentum operator p 9 , can be seen using the commutator identity 



In the context of the matrix element we have the standard result 

_ ifx(W f - Wi) 



(% I q | V4> , 



(40) 

(41) 
(42) 



where Wi and Wf are hydrogenic-type eigenenergies of the initial and final hydrogenic-type states participating in 
the transition, such that h~ 1 (Wf — Wi) = u>. It is clear that the interaction in the dipole approximation involves only 
the center of mass cylindrical coordinates {pr, 4>r, z r)- Substituting this into the transition matrix element we can 
write 



i(Wi - W f ) 



(^ q \e-d\%)U f R (-R);n f i(R) ^(R);n' 



(43) 



where we have replaced q by d = eq, the electric dipole moment vector, and A(R) is the scalar operator of the 
vortex vector potential, the optical polarisation vector e being incorporated into the dipole matrix element, (e ■ d)^. 
Evaluation of the integral of this scalar vortex potential over the initial and final states yields 



M 



fi _ ij4Wf - w;) 



ov 



{¥ q | « • d I V'g) ■A.8 {LtLI _ l) 8 {n i tnt+1) 8{K' z - K z - kg) 



if-i) 



5{K' Z -K z + k z 



(44) 



where the factors A and B arise from the integration, and are not the same. The conservation of OAM is indicated by 
the presence of Kronecker deltas 5(l,l'-i) an d 8(l,L'+1) such that there is an exchange of orbital angular momentum of 
magnitude lh. The first term indicates the process of absorption of a vortex photon by the center of mass, increasing 
the orbital angular momentum of the center of mass about the z-axis by lh, and decreasing the occupancy of the 
photon field by one. The second term describes the reverse process, that of emission of a vortex photon by the 
center of mass. Eq. (|44j) embodies conservation of linear as well as angular momentum. The Dirac delta functions 
6(K' Z — K z — k z ) and 5{K' Z — K z + k z ) exhibit the exchange of linear momentum between the the optical vortex and 
the center of mass such that by absorbing or emitting a vortex photon the linear momentum of the center of mass 
changes by hi k z . 

Thus we conclude that it is possible to transfer orbital angular momentum only between the optical vortex and 
the center of mass, as previously shown by (l8l. [20|. The internal electron- type motion is involved in the process of 
OAM transfer only if the optical vortex beam is circularly polarised, as can be seen by evaluation of the dipole matrix 
element (ip? |e • d\tp q ), but this has nothing to do with the orbital angular momentum of the light field due to the 
vortex factor e %l ^ R . In the dipole approximation, it is not possible to transfer orbital angular momentum from the 
vortex to the internal degrees of freedom of the electron-type motion. This is in agreement with the results of the 
experimental investigations by Araoka et al. In their work, these authors demonstrated that OV light is not 

specific in its interaction with chiral matter, and that orbital angular momentum may only be exchanged between the 
optical vortex and the center of mass of the system. The experimental results also confirm the theoretical predictions 
of earlier investigations using the PZW Hamiltonian approach [20j . Further experimental results have confirmed that 
OAM may be transferred between an OV and the rotational motion of an atom [2~j| . 



IV. THE ELECTRON VORTEX CASE 



For the case of the electron vortex interacting with the two-particle system we have the total Hamiltonian 



H = H° + fl*$, 



(45) 
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where here the zero order Hamiltonian contains only those terms relevant to the two-particle system and the electron 
vortex; 



h° = h2 + hL + M 



EV- 



The interaction Hamiltonian H l £y is given by the Coulomb interaction, expanded to dipolc order: 



/Vint 
-^EV 



q • (r„ - R) 



47re |r„ 

as derived above. The relevant matrix element in this case is 



R| ; 



M 



fi 
EV 



^ (q) ; V>£ (R) ; $1 (r„ ) H% ^ (q) ; ^(R) ; V* (r„) ) , 



which can be rewritten as the following product of matrix elements 



M 



fi 
EV 



R 



Rl 3 



^(r„);^(R) 



(46) 



(47) 



(48) 



The first term can be expressed as the dipole matrix element (ipf | d | f/'g)) where, as before, the atomic electron dipolc 
moment is given by d = eq. The second term in Eq. (|48p. when evaluated, yields 



R 



Rl 3 



C 



fi 



>[(L+l),(L'+l< + l)] 



V 



S[(L+l),{L'+l'-l)] +ZzS[(L+l),(L>+t')], 



(49) 



where C, T> and I are functions which do not contain angular variables and the carets indicate Cartesian unit vectors 
(See Appendix B for details). When the dot product in Eq. (|4"8")l is straightforwardly evaluated, the relevant matrix 
elements of the dipole moment components between internal states are those of the form (q x ± iq y )/2 and q z . These 
would result in the familiar selection rules in optical dipole transitions. The overall matrix element has the form 



M^ v = Q5[(L+l),(L'+l> + l)]5[m,m'-l] + <S%L+i),(£'+i'-l)]#[m,m'+l] + U8[(L+l),(L'+l')]8[m,m'] , 



(50) 



The first term in Eq. (|50[) indicates the possibility of a single unit of orbital angular momentum being transferred from 
cither the center of mass or the electron vortex to the atomic electron, decreasing the total OAM of the combined 
vortex-center of mass system by one unit, and increasing the magnetic quantum number of the atomic electron by one 
unit. The second term describes the reverse process - that of a unit of OAM being transferred from the internal atomic 
motion to the vortex-center of mass system; the third term indicates the possibility of an interaction in which no 
OAM is transferred between the three sub-systems. This is clearly different to the case involving the optical vortex, in 
which only transfer between the vortex and the center of mass is possible, unless the beam itself is circularly polarised. 
This result was previously reported in [22| . 



V. DICHROISM IN EELS 



Predictions of the OAM selection rules for the EV case can be readily put to the test by consideration in relation to a 
very recent experiment by Verbeeck et al. [|| , who set out to measure the dichroic signal in the case of L2 and L3 edges 
in magnetized iron thin film. Their experiment essentially involves a single particle transition matrix element, so only 
the m-selection rule is of importance and it is easy to see that the same OAM selection rules would apply if the single 
particle internal atomic wavefunction in the model system discussed here is replaced by the many-particle electron 
wavefunctions of the transition metal atom (23[. This is accomplished by replacing |^*) = \n\ l\ m) by \2p 6 3d n ; j; ntj), 
and \4>q) = \n;£';m') by \2p 5 3d n+1 ;j';m'j), where j and rrij are the total angular momentum quantum number and 
associated magnetic quantum number of the many-particle transition metal atom in which the 2p core electrons and 
3d valence electron states are involved in the transition [24[ . 

By Fermi's golden rule, the azimuthal angular dependence of the transition rate T±i for the transition process 
involving by an EV mode having I = ±1, is such that 

T ±t oc IM^fpf, (51) 

where pf is the density of final states. As before, we need only consider the angular dependence, since the primary 
purpose is to examine the role of the OAM of the EV mode in the transition process. From Eq. (f5T)|)tfie internal atomic 



9 



dipolc transitions in question are those for which to'- = rrij ± 1 , for application to the result obtained by Verbeeck et 
al. [H]. These transitions can occur accompanied by the transfer of one unit of OAM of I = ±1, either gained or lost 
by the EV mode. Since this interaction takes place within iron thin films, a first approximation is to assume that 
the atomic centers are fixed, ensuring L = L' = 0, so that all OAM exchange occurs only between the vortex and the 
internal motion. The dichroism signal is proportional to the difference between the transition rates for the scenario 
involving the vortex with OAM equal to I on one hand and that involving the vortex with OAM equal to — I on the 
other. Relaxing the fixed atom approximation would involve center of mass states in the form of phonons endowed 
with OAM, which would then participate in the OAM exchange. 

By examining the possible transitions in the L2 and L3 edges we can determine whether or not we expect to see 
dichroism based on the transitions the I = ±1 beams may excite (without, at this point, taking into account density 
of available final states). We deal here with many-electron wavefunctions, and the magnetic quantum number affected 
here is the total angular momentum magnetic quantum number, mj. The iron L edge corresponds to transitions 
from the 2p states to 3d states. This satisfies the selection rule Al = +1, and interaction with the beam requires 
Arrij = ±1. Since we deal with the dipole approximation, we have the restriction that As = 0. Using a very simple 
model in which the angular part of the many-electron wavefunction is expressed in terms of spherical harmonics of 
the form Y™ 3 [25[, we have now a set of 12 allowed transitions - six each allowed for the different beam orbital 
angular momentum polarisations, with these six further subdivided into two from the L2 edge and four from the L3 
edge for both senses of rotation of the beam (see Appendix [B| . Each possible transition in the interaction of the Fe 
atom with the I — +1 beam has a corresponding transition induced by the I = — 1 beam, and the strengths of these 
two interactions are the same (full details are given in Appendix [Bjl . The total transition rate r±i for the L edge is 
given by the sum of the matrix elements over the set of possible intial and final states, from L2 and L3, multiplied 
by the densities of the final states. Thus, since each I = +1 transition has a corresponding I = — 1 transition, and 
the strengths are the same, the square modulus of the matrix elements will be the same for interactions with EVs 
with I = ±1. The dichroic signals will arise from the difference in the to- level occupancy. It should be noted that 
here the m-level occupancy does not refer to that of the 3d levels of the transition metals, but the multiplct states of 
the 2p 5 3d n+1 . Because of the strong spin-orbit coupling between the core hole and the rest of the electronic systems, 
these signals will be sensitive to the spin magnetization as well 24 1 . 

Verbeeck et al. measured a clear dichroism signal in their experiment Q. Using the analysis above, this can be 
explained as due to the difference in the density of states such that pf ^ pj. Since these processes are dominated by 
dipole active transitions the dichroism signal can be detected using small angle scattering. Clearly this is advanta- 
geous when compared with the electron energy loss magnetic chiral dichroism technique used by Schattschncider et 
al. pr| . which operates using large scattering angles, and where the cross-section for energy loss is much reduced and 
experiments arc more susceptible to noise. 



VI. COMMENTS AND CONCLUSIONS 



We have demonstrated by explicit analysis that a transfer of OAM can indeed occur between an EV mode and the 
internal electronic-type dynamics of matter involving an electric dipole transition. This contrasts sharply with the 
case of optical OAM exchange in the interaction of an OV mode with similar systems. Our predictions that optical 
vortices cannot exchange OAM with the internal dynamics of an atom in the leading electric dipole interaction are 
consistent with the experimental finding of Araoka et al. (l6j who demonstrated that OVs are not specific in their 
interaction with chiral molecules. This is also consistent with earlier theoretical predictions by Babiker et al. [2(| who 
used a complicated PZW technique to analyse the interaction. This significant finding implies that the EV beam 
should play a more effective role in magnetic energy loss spectroscopy than that played by ordinary electron beams 
and, equally significantly, than the role played by light beams. It is well recognised that the dipole transition is 
often the dominant process in most physical systems. In particular, it should now be possible for EV beams to be 
used to detect circular dichroic activity in p roteins and other biological molecules, allowing useful information to be 
gained about their secondary structures [27l l28j . It is expected that the sensitivity and spatial resolution will be high, 
provided that radiation damage effects can be mitigated. 

In view of the similarities of EVs and OVs it is clearly natural to contemplate whether an experiment similar to 
that by Araoka et al. [l6| on the handedness of processes involving the possibility of OAM exchange with internal 
dynamics should be carried out for the EV case. As far as the authors are aware, the only experiment to date involving 
the transfer of OAM of electron vortices is that by Verbeeck et al. Q who specifically dealt with the case I = ±1 to 
investigate the electron energy loss spectroscopy signals from magnetised Fe films using EV beams. The transitions 
involved are those by core electrons participating in electric dipole allowed transitions between discrete atomic energy 
levels within a rigid condensed matter background of essentially fixed Fe atoms. The expectation was that the two 
EELS signals, one from I = 1 and the other from / = —1, would be different, suggesting that the EELS revealed 
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an intrinsic chirality of the medium. This was indeed the case in the experiment by Verbeeck et al. The question, 
however, arises as to the theoretical basis for the observed dichroism. 

Our task in this paper has been two-fold. Firstly, we set out to construct the basic theory underlying the coupling 
of an OV and an EV to matter represented by a two-particle atom as a neutral bound system of two charged particles. 
The basic theory has been carried out in a manner which incorporates both types of vortex and as far as the atom is 
concerned, it was necessary to distinguish internal (electronic-type) and gross (center of mass-type) degrees of freedom 
of the atomic system. Secondly we proceeded to consider each vortex separately in its interaction with the atomic 
system. We have shown that, contrary to the case of OVs, the theory for EVs allows the transfer of OAM to the 
internal dynamics for electric dipolc transitions, as in the experiment by Verbeeck et al. It turns out that this crucial 
difference between EVs and OVs is attributed to the distinct interaction mechanisms of the vortex with matter. 

Our results suggest that the matrix elements involving the orbital angular momentum transfer has the general 
vector product form of (eq)^ • (f t ,(r))y,- where (eq)^ is the dipole transition matrix element between internal atomic 
states and f«(r) is the effective field seen by both the internal motion and center of mass of the atomic system. It 
turns out that the condition for the transfer of orbital angular momentum between the vortex beam and the internal 
motion of the atom is for f v (r) to exhibit chirality. In the case of EVs, the longitudinal Coulomb interaction couples 
the dipole moment of the internal motion of the atom to the electric field of the vortex beam, in a similar manner to 
that in which atomic electrons couple to the optical field of circularly polarised light. In the case of OVs, in which 
the vortex is characterised by a linearly polarised transverse vector potential carrying orbital angular momentum, the 
transition matrix element only depends on the value of the transverse vector potential at the atomic site, hence the 
dot product is not chiral-spccific and no OAM transfer to the internal dynamics is allowed to accompany an electric 
dipole transition. A transfer of OAM to both parts of the atomic system in higher multipolar transitions than the 
dipolar is, however possible for both types of vortex. 

Additionally, we note that the transfer of OAM from the EV to the atom may proceed independent of the dynamics 
of the atom, due to the long range Coulomb interaction involved. This is in contrast to the case involving OV beams 
that requires the center of mass of the atom to be a dynamical variable, free to rotate about the beam axis [l8| . 

We have analysed the theory for the general case of topological charge I. However for the specific case of the 
Verbeeck et al. experiment, we have demonstrated that although OAM exchange can occur between the EV and 
matter in electric dipole transitions for the opposite helicitics I = ±1, there is no intrinsic difference in the EELS 
absorption between the two helicitics. We have shown that the experiment by Verbeeck et al. displays a dichroism 
effect via a similar mechanism to the analogous XMCD experiment [24| . 
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Appendix A: EV matrix element 



Starting from 
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' V 
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^(r 9 ) 



i>l(r q ) • V, J (r,);^( R ) 



R 



R 



^(r,);Vk(R) 



(El 



we seek to fully evaluate this matrix element, by considering it as the product of two matrix elements - the first 
being the dipole matrix element between hydrogenic wavefunction, and the second that of the Coulomb potential 
between the center of mass states and vortex wave functions. The dipole matrix element is well known, so we will first 
concentrate on this second, Coulombic matrix element. This is best evaluated using Cartesian coordinates, so as to 
compare the result with the known dipole result, and compute their scalar product. Thus, 



(r„ - R) (p v cos(^) - p R cos(<j) R ))}L+ (p„sin(0„) - p R sin((j> R ))y + (z v - z R )i 



|r„-R| 3 



{F{p v , z v , p R , z R ) - G(pv,p R ) cos((j) v - 4> R )Y 



(Al) 



Evaluating this as the matrix element between (r v ); ip R ^(R)) requires the use of the substitution 
in order to express the intergral in Eq. (|48[) in terms of generic integrals of the form 



2tt 



J(l-l'+a)y 



-dy, 



(A2) 
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where a is an integer taking the values 0, ±1. The result, after making this substitution, may be written as 

rJ-RI 3 ) = C \2J 6 K L + l )^'+i'+i)] + v (j-2) <5 [( i + ; )^(i'+''-i)] + 22 < y [(i+0.(i'+i')] ® 
where 

c = K y. 1 -xy , (A3) 

v = K y +1 -xy , (A4) 

I=vyo, (A5) 

with k, A and 77 being factors arising from the integration over the remaining (non-azimuthal) degrees of freedom of 
ip v (r v ) and i/jr(R). 

Combining the two parts involved in the dot product of the overall transition matrix elelment we have 



fi 

V 

+ y (?x - iq V ) fi %L+Z),(L'+Z'-1)] + 1 fi <%L+Z).(L'+f')] ( A6 ) 

The evaluation of each of the dipole component matrix elements ^ (q x ± iq y ) ^ and (q z )if is standard and leads to the 
usual m-selection rules. The overall matrix element can be written as 

-M E y = Q5[(L+l),(L'+l' + l)]3[ m ,7n'-l] +<S%L+0,(-£/+i'-l)]£[m,m'+l] + U5[(L+l),(L'+l')]$[m,m'] , fSD]) 

where Q, S and W are factors containing the normalisation constants and integrals not involving azimuthal angles. 

Appendix B: Application of the selection rules 

So as to compare our results with the experimental results of we now apply the selection rules derived for the 
electron vortex interaction with the two-particle system, by considering transitions between certain initial and final 
states. The transitions involved are those of core electrons participating in electric dipole allowed transitions between 
discrete atomic energy levels within a rigid condensed matter background of essentially fixed Fe atoms at equilibrium. 
If the center of mass of the atoms in the lattice is allowed to participate in the angular momentum transfer, the 
motion would involve phonon states endowed with OAM. The electron vortex interaction matrix element is 

M^y = QS[(L+t),(L'+l'+l)]fi[m,m'-i\ + £%L+Z),(L'+Z'-l)]<5[m,m'+l] + ^%L+Z),(I/+Z')]^[m,-m']j CUB 

The factors Q and S can be written as follows 

Q = K'C = K' (k3>_i - \y ) 

s = kv = k { K y +1 - xy ) . 

where K. and KJ arise in the evaluation of the dipole matrix elements 5 (q x — iq y ) i j. From the definition of y a we can 
compare these intensity factors for a specific transition. Choosing I = ±1, and I' = we can write 

/ /-27T gO /.27T i(„) \ 

C=U/ sdy-X -dy), (Bl) 

\ Jo (T-Gcos(y)) 2 Jo (J 7 - G cos(y)) 2 / 

2tt e ^2tt e -i(y) 



V = \ k -dy - X / -dy , (B2) 

\ Jo (T-Gcos(y)) 2 Jo {T-Gcos(y)) 2 J 

From which it can be seen that C = V* . This choice of I and V suggests m' = m+ 1 for I = 1 and m' = to — 1 for I = — 1 , 
from some initial magnetic state to. This information can be used to find the factors /C, K! for each case. This now 
needs generalising to the many-electron wavefunction of the model iron atom. In the LS coupling regime, the angular 
part of the wavefunction is given by the product of the spherical harmonics of the occupied states. The total orbital 
angular momentum of the atom is J = L + S, and associated magnetic quantum number, rrij, is the quantity that 
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TABLE I. Allowed transitions in the L2 and L3 edges for interaction with an electron vortex beam with / = +1 and I = — 1. 



will be affected by the absorption or emission of a unit of OAM from the electron vortex. We can express the angular 
part of the many-electron wavefunction in terms of spherical harmonics of the form Y™" 3 (neglecting numerical and 
phase factors that arise from the coupling of the constituent electrons) [25| . The spherical harmonics are normalised 
such that 

(-l) e -^Yf m = Yf 1 *, (B3) 

and this also applies to our many-electron spherical harmonics Yj j . For our purposes, we seek to explain the dichroism 
of the L-edge transitions of magnetised iron thin films observed in EELS using electron vortex beams, by Verbeeck 
ct al. To do this, we first look at the possible transitions that can be excited in the model iron atom by the 
two beams of opposite orbital angular momenta, I = ±1. These are summarised in Table H] As can be seen, each 
transition in the / = +1 case has a corresponding, similar, transition in the I = —1 case, having rrij = — m^- and 

m'j 1 ' = —m'j 1 '. It is clear that, due to Eq. (|B3[) . the dipole matrix elements of the corresponding transitions will 
have the same magnitude (due to the symmetry properties of the Wigner 3-j symbols used to calculate the coupled 
spherical harmonics, the strengths for the forward and reverse transitions for I = ±1 remain the same (23[, so the 
numerical factors can be safely neglected, as stated above). 

The total transition rate (as observed in Q ) of the L 2 and L 3 edges will be given by the sum of the transition rates 
of the individual transitions in each edge. So, for the L2 edge we have 



|C| 2 (j <2pi/2,Tij = -l/2|q|3d 3/ 2,mj = +1/2) \ 2 p 3 d 3/2 ( mj =+i/2) 

+ I (2p 1/2 ( mj = +1/2) I q I 3^3/2,77^ = +3/2) | 2 /53d 3 / 2 (™ J =+3/2)^) , (B4) 



and 



It-" 1 = Y^ 2 i ( 2 P^> m l = +1/2 I q| 34/2,77^- = -1/2) r / 33d 3/2 (m 3 .=-l/2) 



+ I (2pi/2, rrij = -1/2 | q | 3d 3/2 , rrij = -3/2) | 2 P3d 3/2 ( mj =~3/2)^ , (B5) 

which will be equal as long as p3d 3/2 ( mj =+i/2) = P3d a/2 ( mj =-i/2) and p 3 d 3/2 , mj =+3/2 = P3d 3f2 ., m] =-3/2, as we have 
already established that \C\ 2 = \T>\ 2 . The same argument applies to the L3 edge. Thus we conclude that the observed 
dichroism is due solely to the distribution of electrons in the magnetised iron, and not due to the mechanism of the 
interaction with the electron vortex. 
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